AD-A150  885  PRESSURE  DISTRIBUTION  ON  PROPELLER  BLADE  SURFACE  USING  i/i 
NUMERICAL  LIFTING  (U)  DAVID  U  TAVLOR  NAVAL  SHIP 
RESEARCH  AND  DEVELOPMENT  CENTER  BET  K  H  KIM  ET  AL. 
UNCLASSIFIED  JAN  85  DTNSRDC-84/072  F/G  20/4  NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS-1963-A 


HFPnonucrn  at  government  expense 


•f;--'  *V„ 


DAVID  W.  TAYLOR  NAVAL  SHIP 
RESEARCH  AND  DEVELOPMENT  CEI 

,  i.  ./  .  • 


Bathaada,  Maryland  20084-5000 


•  -s.  . 

•> 

;  V"  .  ; 

PRESSURE  DISTRIBUTION  ON  PROPELLER  BLADE  SUR^E  ' 

USING  NUMERICAL  LIFTING  SURFACE  THEORY  ■ 

••  • 


...'V.* /..>v 

' 


Ki-Han  Kim 
Sukeyuki  Kobayashi 


APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UlftlMITED 


Presented  at  Propellers  1984  Symposium  V  '  ' 

Virginia  Beech,  Va.  . 

15-16  May  1984 

■  w 


SHIP  PERFORMANCE  DEPARTMENT 
RESEARCH  AND  DEVELOPMENT  REPORtr^^^g 


January  1908tl.A~'^R^**'-- 


• 

r 

/  M 

V;  .  •  t,!-" 

» 

1  • 

r-_ 

'  • 

» 

.  _  .  . 

.  ... 

UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  RAGE  (When  Bata  Entered) 

REPORT  DOCUMENTATION  PAGE 

r  RERORT  NUMBER  Jj.  GOVT 

DTNSRDC-84/ 072 

«.  TITLE  (end  Submit)  9 


J2.  GOVT  JwC  CESSION  NO. 


READ  INSTRUCTIONS 

_ BEFORE^gQMPLETINC  FORM 

s^RECirirNt's  catalog  number 


S.  TYPE  OF  REPORT  S  PERIOO  COVEREO 


PRESSURE  DISTRIBUTION  ON  PROPELLER  BLADE 
SURFACE  USING  NUMERICAL  LIFTING 
SURFACE  THEORY 
~  autmorc*;  ~ 

Ki-Han  Kim 
Sukeyuki  Kobayashi 

*  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

David  W.  Taylor  Naval  Ship  Research 
and  Development  Center 
Bethesda,  Maryland  20084-5000 

"  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

David  W.  Taylor  Naval  Ship  Research 
and  Development  Center 
Bethesda,  Maryland  20084-5000 

I  '«  MONITORING  AGENCY  NAME  A  ADDRESS?/?  different  Iron,  Controlling  Otttce) 


I  16  DISTRIBUTION  STATEMENT  (ot  thle  Report) 


Final 

S.  PERFORMING  ORG.  REPORT  NUMBER 


•  •  CONTRACT  or  GRANT  NUMBERO) 


10.  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 


12.  REPORT  DATE 

January  1985 

13.  NUMBER  OF  PAGES 

54 

Vs.  SECURITY  CLASS,  (ol  thle  report) 


UNCLASSIFIED 


ISa.  DECLASSIFICATION' DOWNGRADING 
SCHEDULE 


APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UNLIMITED 


[  17.  DISTRIBUTION  STATEMENT  (ol  the  eb.tr  ect  entered  In  Block  20,  II  dlllerent  from  Report) 


10  SUPPLEMENTARY  NOTES 


Presented  at  Propellers  1984  Symposium, 
Virginia  Beach,  Virginia,  15-16  May  1984 


'»  KEY  WORDS  (Continue  on  reveree  aid#  II  neceeeery  and  Identity  by  block  number)  " 

Propeller.  ‘  4  —i  )  ■  L 

cBTade  Pressure./ 

Lif  ting  „  Surf  ace  Theory > 

Discrete  Vortex/Source  Method  - 

20.  ABSTRACT  (Continue  on  reverie  elde  II  neeeeeary  and  Identity  by  block  numbet)  ’ 

yA  procedure  and  numerical  results  are  presented  for  the  prediction  of 
the  steady  pressure  distribution  on  a  rotating  propeller  blade  surface 
based  on  lifting  surface  theory.  A  computer  code,  named  the  Propeller 
Steady  Pressure  (PSP)  program,  has  been  developed  by  extending  the  existing 
propeller  analysis  program,  PSF,  based  on  vortex/source  lattice  techniques, 
developed  at  the  Massachusetts  Institute^  of  Technology.  Predictions  by 

(Continued  on  reverse  side) 


00  1473? 


EDITION  OF  1  NOV  ••  It  OBSOLETE 
$/N  0102-  LF-  014-  6601 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  RAGE  i 


Ban  Data  Bntered) 


UNCLASSIFIED _ 

SECURITY  CLASS!  TIC  ATIOR  Of  THIS  RARE  (9km  Ml  fclwi 

(Block  20  continued)  — , 

-PSP  are  compared  with  selected  experimental  values  that  are  believed  to 
be  accurately  and  reliably  Measured.  Comparisons  are  also  made  between 
PSP  predictions  and  other  theoretical  predictions.  The  predictions  by 
PSP  are  generally  in  good  agreement  with  experimental  values  and  with 
other  prediction  methods  except  for  the  tip  region  where  current  pro¬ 
cedures  may  not  be  accurate  enough  to  represent  the  actual  flow. 


TABLE  OF  CONTENTS 


Page 


LIST  OF  FIGURES .  iii 

TABLE .  iv 

NOTATION  .  v 

ABSTRACT  .  1 

INTRODUCTION  .  1 

OVERVIEW  OF  PSP . 2 

MODIFICATION  TO  PSF .  4 

RESULTS  AND  DISCUSSIONS  .  8 

TWO-DIMENSIONAL  AIRFOIL  SECTIONS  .  8 

IHI  MODEL  PROPELLER  MP  282  .  8 

DTNSRDC  PROPELLER  4718  .  10 

NSMB  MODEL  PROPELLER .  11 

COMPARISON  WITH  OTHER  PROCEDURES  .  11 

SUMMARY  AND  CONCLUSIONS .  12 

RECOMMENDATIONS  .  13 

ACKNOWLEDGMENT .  13 

REFERENCES .  15 

APPENDIX  -  PRESSURE  IN  A  MOVING  FRAME  OF  REFERENCE  .  17 


LIST  OF  FIGURES 

1  -  Fixed  and  Rotating  Coordinate  Systems  for  a  Right-Hand- 


Rotation  Propeller .  21 

2  -  Discretization  of  Blade  Singularities  .  22 

3  -  Horseshoe  Vortices  on  Blade  and  in  Wake .  23 

4  -  Trailing  Wake  Geometry  after  Wake  Alignment  .  24 

5  -  Discrete  Singularity  Distribution  for  Two-Dimensional 

Airfoil  Section  .  25 


6  -  Schematic  Representation  of  the  Effect  of  Chordwise  Vortices 

on  the  ich  Pressure  Point  . .  26 

7  -  Pressure  Distribution  on  Flat  Plate  at  a  =  4  Degrees .  27 

8  -  Pressure  Distribution  on  NACA  0012  Section  at  a  =  0 

Degrees . 28 

9  -  Pressure  Distribution  on  NACA  4412  Section  at  a  =  6.4 

Degrees .  29 

10  -  Open-Water  Performance  of  IHI  Model  Propeller  MP  282  .  30 

11  -  Pressure  Distribution  on  IHI  Model  Propeller  MP  282 

at  J  =  1.054  31 

12  -  Pressure  Distribution  on  IHI  Model  Propeller  MP  282 

at  J  =  1.163 . 32 

13  -  Surface  Flow  Patterns  by  Oil-Film  Test  on  IHI  Model 

Propeller  MP  282  33 

14  -  Pressure  Distribution  on  DTNSRDC  Propeller  4718  at 

J  =  0.75 .  34 

15  -  Pressure  Distribution  on  Suction  Side  of  DTNSRDC  Propeller 

4718  at  r/R  =  0.5:  Comparison  of  Prediction  and  Two 

Sets  of  Experiments .  35 

16  -  Open-Water  Performance  of  NSMB  Model  Propeller  .  36 

17  -  Pressure  Distribution  on  NSMB  Propeller  at  J  =  0.4 .  37 

18  -  Pressure  Distribution  on  NSMB  Propeller  at  J  =  0.6 .  39 

19  -  Pressure  Distribution  on  DTNSRDC  Propeller  4118  at 

J  =  0.833:  Comparison  with  ITTC  Propeller 

Committee  Survey  .  41 

20  -  Pressure  Distribution  on  DTNSRDC  Propeller  4498  at 

J  =  0.888  .  42 


Table  1  -  Propeller  Geometric  Characteristics . . .  43 


iv 


NOTATION 
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J  =  V/nD 

Kq  =  Q/(pn2D5) 
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n 
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Blade-section  chord  length 

Pressure  coefficient 

Propeller  Diameter 
Meanline  shape  function 
Gravitational  acceleration 

Unit  base  vectors  in  a  cylindrical  polar  reference  system 

Unit  base  vectors  in  a  Cartesian  reference  frame 

Total  rake:  axial  displacement  of  blade-section  midchord 
point  from  y-z  plane 

Advance  coefficient 
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Thrust  coefficient 

Propeller  votational  speed,  revolutions  per  unit  time 

Pitch  of  blade  section 

Pressure 
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ABSTRACT 


A  procedure  and  numerical  results  are  presented  for  the 
prediction  of  the  steady  pressure  distribution  on  a  rotating 
propeller  blade  surface  based  on  lifting  surface  theory.  A 
computer  code,  named  the  Propeller  Steady  Pressure  (PSP) 
program,  has  been  developed  by  extending  the  existing  pro¬ 
peller  analysis  program,  PSF,  based  on  vortex/source  lattice 
techniques,  developed  at  the  Massachusetts  Institute  of 
Technology.  Predictions  by  PSP  are  compared  with  selected 
experimental  values  that  are  believed  to  be  accurately  and 
reliably  measured.  Comparisons  are  also  made  between  PSP 
predictions  and  other  theoretical  predictions.  The  pre¬ 
dictions  by  PSP  are  generally  in  good  agreement  with  experi¬ 
mental  values  and  with  other  prediction  methods  except  for 
the  tip  region  where  current  procedures  may  not  be  accurate 
enough  to  represent  the  actual  flow. 


INTRODUCTION 


Knowledge  of  the  pressure  distribution  on  the  propeller  blade  surface  is 


essential  to  understanding  cavitation  phenomena,  boundary  layer  characteristics 
and  stress  on  blades.  Measuring  the  pressure  distribution  on  a  rotating  blade  is 
extremely  difficult  and  time-consuming,  and  even  then  the  reliability  and  repeat¬ 
ability  of  the  experimental  data  are  often  questionable.  Nevertheless,  a  number 


of  experimental  results  of  reasonable  reliability  are  available,  such  as  the  ones 


1*  2  3  4  5 

obtained  by  Mavludoff,  Kato,  Yamasaki,  Takei  et  al.,  Jessup,  and  Versmissen 


and  Van  Gent. 


The  ability  to  predict  the  blade  pressure  distribution  reliably  and  accurately 


is  also  highly  desirable.  Many  institutions  throughout  the  world  have  their  own 


prediction  methods;  most  of  them  are  based  on  lifting-surface  methods  such  as  those 


7  8  9  10 

of  Okamura,  Kuiper,  Brockett,  and  Tsakonas  et  al.  and  a  two-dimensional  pro¬ 


cedure  with  some  empirical  corrections  for  three-dimensional  effects  by  Bahgat. 

In  this  report,  a  procedure  is  presented  to  predict  the  pressure  distribution 


on  the  propeller  blade  surface  operating  in  steady  flow  based  on  the  discrete 

12 

vortex/source  lattice  method  developed  by  Greeley  and  Kerwin.  The  discrete 


vortex  lattice  method  has  been  used  in  the  field  of  aerodynamics  as  early  as  1943 


*A  complete  listing  of  references  is  given  on  page  15. 
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by  Faulkner  for  the  calculation  of  aerodynamic  forces  on  an  arbitrary  wing  shape. 

14 

The  accuracy  of  this  simple  method  has  been  found  very  satisfactory  and  in  two- 

dimensional  flow  "remarkable”.^  The  primary  advantage  in  using  the  discrete 

vortex/source  lattice  method  is  the  ease  and  the  flexibility  to  model  the  complex 

geometries  of  the  propeller  blades  and  their  trailing  vortex  wake.  With  the  advent 

of  large  computers,  panel  methods  are  widely  used  for  the  design  and  analysis  of 

three-dimensional  aerodynamic  configurations  both  as  the  simple  vortex/source 

lattice  approximation  and  as  more  complex  local  elements. 

In  the  area  of  marine  hydrodynamics,  Kerwin  and  Lee^  developed  a  discrete 

vortex/source  method  and  corresponding  computer  code,  PUF2,  for  the  prediction  of 

steady  and  unsteady  performance  of  subcavitat ing  propellers.  Kobayashi^  and 
1 8 

Kim  developed  a  procedure  to  compute  the  pressure  distribution  based  on  the 

16  12 
method  developed  by  Kerwin  and  Lee.  More  recently,  Greeley  and  Kerwin 

developed  design  and  analysis  procedures  and  corresponding  computer  codes,  PBD-10 

for  design  and  PSF  for  analysis,  for  propellers  operating  in  steady  flow.  Greeley 
12 

and  Kerwin  made  improvements  to  the  steady  part  of  the  procedure  developed  by 

Kerwin  and  Lee*^  in  two  major  areas;  one  is  the  improved  semi-empirical  description 

of  the  trailing  vortex  sheet  and  the  other  is  the  capability  to  model  the  flow 

over  the  outer  portion  of  the  blade  more  accurately.  In  the  present  work,  only 

1 2 

the  "global"  part  of  the  procedure  developed  by  Greeley  and  Kerwin*  has  been 
investigated  for  the  pressure  distribution.  The  more  accurate  local  flow  model 
near  the  tip  is  yet  to  be  examined. 

A  computer  code,  named  the  Propeller  Steady  Pressure  (PSP)  program,  has  been 

1  9 

developed  by  extending  the  propeller  analysis  program,  PSF,  presented  by  Greeley 

12 

and  Greeley  and  Kerwin.  This  report  describes  the  computer  code  and  presents 
some  comparisons  of  the  predictions  made  by  PSP  with  experimental  measurements  and 
predictions  by  other  theories. 


OVERVIEW  OF  PSP 

The  Propeller  Steady  Pressure  (PSP)  code  is  basically  the  same  as  the  Propeller 

12 

Steady  Flow  (PSF)  analysis  program  developed  by  Greeley  and  Kerwin,  except  for 
the  additional  capability  in  PSP  of  calculating  the  pressure  distribution  on  the 
blade  surface.  The  PSF  code  assumes  that  the  propeller  operates  in  an  axisymmetric 
onset  flow  consisting  of  axial,  radial,  and  tangential  components.  The  presence 
of  the  propeller  hub  and  any  other  boundaries  to  the  flow  is  ignored.  The  blade 


2 


boundary  layers  are  assumed  to  be  thin,  so  that  the  flow  can  be  treated  as  inviscid 
except  for  the  calculation  of  frictional  drag. 

The  nonrotating  coordinate  system,  (xo,yo,zQ),  and  rotating  coordinate  system, 
(x,y,z),  fixed  to  the  blades  are  shown  in  Figure  1.  The  x-axis  of  the  fixed  and 
rotating  system  are  coincident,  as  are  the  (y,z)  and  (y  ,z  )  planes.  The  defini¬ 
tion  of  the  angular  coordinates  in  the  fixed  system,  0  ,  and  in  the  rotating  sys¬ 
tem,  0,  are  also  defined  in  Figure  1.  The  propeller  rotates  at  a  constant  angular 
velocity,  ft  =  -£H.  A  field  point,  P,  in  the  fluid  with  angular  coordinate,  0,  in 
the  rotating  frame  has  an  angular  coordinate 

0  =  6  -  ftt  (1) 

o 

in  the  fixed  frame  for  a  right-handed  propeller  shown  in  Figure  1. 

The  blade  geometry  is  defined  relative  to  a  midchord  line,  which  is  para¬ 
metrically  defined  by  the  radial  distribution  of  skew,  0g(r),  and  total  rake, 
i,„(r).  The  pitch  angle,  4>  ( r ) ,  and  chord  length,  c(r),  define  the  angle  and  extent 
of  the  sectional  nose-tail  line  along  the  pitch  helix  on  the  surface  of  a  cylinder 

of  radius  r.  The  meanline  offset,  f(r,x  ),  and  thickness  distribution,  t(r,x  ), 

c  c 

describe  the  section  characteristics  of  the  blade  as  a  function  of  radius,  r,  and 
nond imens iona 1  arc  length,  xc,  along  the  nose-tail  line.  The  meanline,  f,  is 
measured  along  the  cylindrical  surface  at  right  angles  to  the  nose-tail  line.  The 
thickness,  t,  is  measured  perpendicular  to  the  meanline.* 

The  blades  and  vortex  wake  are  represented  by  straight-line  vortex  and  source 
lattice  elements  of  constant  strength,  distributed  over  the  meanline  surface  of  the 
blade  (see  Figure  2)  and  the  assumed  surface  of  the  trailing  vortex  sheet.  The 
vortices  are  arranged  in  the  traditional  horseshoe  configuration  (see  Figure  3)  so 
as  to  satisfy  Kelvin's  conditions  automatically,  and  the  strength  of  each  horseshoe 
vortex  is  determined  by  solving  a  set  of  simultaneous  equations,  each  satisfying 
the  flow  tangencv  condition  at  a  blade  control  point.  Source  strength  is  deter¬ 
mined  from  the  slope  of  the  thickness  distribution  and  resultant  onset  speed. 

*  At  DTNSRDC,  the  thickness  is  conventionally  measured  perpendicular  to  the 
nose-tail  line.  In  linear  theory  the  differences  of  these  two  specifications  is 
of  higher  order. 


The  position  of  the  shed  vortex  sheet  is  determined  iteratively  by  first 
solving  the  boundary  value  problem  with  an  assumed  position,  and  then  aligning  the 
wake  with  the  computed  total  velocity  field  for  a  specified  radial  contraction. 

The  boundary  value  problem  is  then  re-solved  and  the  procedure  is  repeated  until 
convergence  (see  Figure  4).  This  process  of  wake  alignment  is  different  from  the 
simple  wake  model  in  PUF2,^  where  the  trailing  vortex  wake  geometry  is  defined  at 
the  outset  by  several  serai-empirically  determined  geometric  parameters. 

Once  a  converged  solution  is  obtained,  blade  forces  are  computed  by  applying 

the  Kutta-Joukowski  and  Lagally  theorems.  The  Lagally  theorem  is  used  to  compute 

the  forces  on  the  source  elements  as  a  modification  for  the  effect  of  the  thickness 

(source).  This  modification  is  equivalent  to  subtracting  the  thickness-induced 

velocity  from  the  total  velocity  used  to  compute  the  Kutta-Joukowski  force  on  the 

vortex  elements.  If  the  thickness-induced  velocity  were  included  in  the  total 

velocity,  the  resulting  Kutta-Joukowski  force  would  be  larger  than  experimental 

values.  In  PSF,  as  in  PUF2,  an  empirical  suction  factor  is  used  to  estimate  the 

leading-edge  suction  force  at  off-design  conditions.  The  reader  is  referred  to 
12 

Greeley  and  Kerwin  for  details  of  the  computation. 

MODIFICATIONS  TO  PSF 

In  PSF,  the  overall  blade  load  is  computed  by  summing  up  the  elementary  loads 
(the  jump  in  pressure  across  the  surface)  acting  on  each  line  vortex  and  source 
element.  The  elementary  load  is  computed  at  the  midpoint  of  each  spanwise  and 
chordwise  singularity  on  the  key  blade  by  assuming  the  average  velocity  over  the 
length  of  a  singularity  can  he  approximated  by  the  velocity  at  its  midpoint.  This 
point  is  called  "load  point."  Since  the  total  velocity  is  calculated  at  each  load 
point  to  compute  the  load,  it  is  logical  to  choose  the  same  point  as  the  "pressure 
point"  for  pressure  calculation.  In  the  present  study,  pressure  is  computed  at 
only  the  pressure  points  on  the  spanwise  singularities  and  is  interpolated  at 
specified  radii. 

The  velocity  calculated  at  the  load  point  in  PSF  is  a  mean  velocity  that  does 
not  include  the  self-induced  velocity  due  to  the  singularity  segment  where  the 
elementary  load  is  calculated.  However,  when  computing  the  pressure,  not  the  jump 
in  pressure,  the  velocity  jump  across  the  singularity  must  he  included. 
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Since  the  vortex/source  sheet  on  the  blade  surface  is  represented  by  "discrete" 
singularity  elements,  each  discrete  element  represents  a  certain  area.  Therefore, 
when  computing  the  velocity  jump  across  the  vortex  sheet,  we  have  to  redistribute 
this  concentrated  vortex/source  over  the  area. 

Consider  a  two-dimensional  airfoil  illustrated  in  Figure  5.  The  discrete 
bound  vortices/sources  are  located  on  the  meanline  at  the  quarter  chord  of  each 
meanline  segment  to  approximate  the  continuous  distribution  of  the  vortex/source 
along  the  meanline.  Suppose  I\  is  the  strength  of  the  bound  vortex  at  the  i*”^1 
segment  whose  length  is  Ac...  Then  the  distributed  vortex  strength,  y.,  over  this 
segment  can  be  approximated  by: 


(2) 


assuming  the  vorticity  is  uniformly  distributed  over  the  segment.  The  velocity 
jump  across  the  vortex  sheet  is  related  to  the  local  vortex  strength,  y^,  as 
follows : 


<\>1  ■ 


and 


(3) 


where  the  plus  sign  represents  the  upper  surface  and  the  minus  sign  the  lower  sur¬ 
face.  In  this  two-dimensional  case,  the  velocity  jump  is  tangent  to  the  surface 
in  the  chordwise  direction. 

Similarly,  the  distributed  source  strength,  ,  over  the  same  segment  will  be: 


(A) 


where  is  the  strength  of  the  discrete  source  element.  The  source  sheet  induces 
a  jump  in  normal  velocity,  that  is  related  to  the  local  source  strength,  q^,  as 
follows : 


(V  )T  =  -j- 


and 


(Vi  - 


(5) 
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where  the  plus  and  minus  signs  represent  the  upper  and  lower  surfaces,  respec¬ 
tively. 

For  three-dimensional  flow  such  as  that  on  propeller  blades,  the  direction  of 
velocity  jump  depends  on  both  the  spanwise  and  chordwise  vortices.  In  this  case, 
both  spanwise  and  chordwise  singularities  have  to  be  properly  accounted  for  when 
computing  the  velocity  jump.  The  following  is  the  algorithm  adopted  in  the  pre¬ 
sent  study. 

th  th 

Suppose  we  want  to  compute  the  velocity  jump  at  the  i  pressure  point  on  i 
spanwise  vortex  element.  The  total  distributed  vortex  at  this  point,  y. ,  is  the 
sum  of  the  spanwise  and  chordwise  distributed  vortices: 

Xi  =  (Xs)i  +  (Xc)i  (6) 

The  spanwise  distributed  vortex,  (^)  ,  is  approximated  by: 

(Ig>i 

(Vi  =  -i?7-  <7> 

i 

where  Ac^  is  the  length  of  the  chordwise  segment  represented  by  the  discrete  span- 

wise  vortex,  (T  ) . .  This  is  analogous  to  the  two-dimensional  distributed  bound 
~ s  1 

vortex  (see  Equation  (2)).  The  chordwise  distributed  vortex,  (yO ,  is  approximated 

by  the  vector  average  of  the  four  adjacent  chordwise  vortices,  (T  ) . ,  (T  )„,  (T  )_, 

— c  I  — c  l  — c.  J 

and  (T  ).  (see  Figure  6): 

— c  4 

<lc>i  ■  T  T, 

i=  1 

where  Arn  is  the  length  of  the  radial  segment  represented  by  each  discrete  chord- 
wise  vortex,  (T  )  .  The  total  distributed  vortex,  y.,  is  then  converted  to  the 
velocity  jump  in  the  tangential  direction  by  using  Equation  (3). 

The  velocity  jump  due  to  the  source  sheet  is  identical  to  the  two-dimensional 
case  (see  Equation  (5))  since  the  boundary  condition  for  thickness  effects  results 
in  the  same  relation  between  source  strength  and  slope  of  the  chordwise  thickness 
distribution  with  radius  as  a  parameter.  These  velocity  jumps  due  to  vortices  and 
sources  are  added  to  the  velocity  induced  by  all  other  singularities  to  obtain  the 
total  velocity  induced  by  the  propeller. 


(£) 

— c  n 

Ar 


(8) 
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The  pressure  on  a  propeller  blade  surface  rotating  at  a  constant  angular 
velocity,  ft  =  -fti  (see  Figure  1),  in  an  axisymmetric  onset  flow  can  be  expressed 
as  (see  Appendix) : 

P  =  -  yP{X'^  -  (Vx)j[  -  rV}  +  pA  (9) 

total  velocity;  V  =  V  +  rQe„  +  V  +  V 
7  *  —  -w  ^6  — p  — o 

axisymmetric  onset  flow;  V  =Vi+Ve  +  V_en 
7  — w  x—  r^r  0-6 

perturbation  velocity  due  to  the  propeller  blades  and 
their  wakes 

perturbation  velocity  due  to  the  other  sources  such 
as  appendages  or  lifting  surfaces 

unit  vectors  in  the  axial,  radial  and  tangential 
directions  in  the  cylindrical  coordinate  system 
(x,r,8)  rotating  with  the  propeller 

in  Equation  (9)  indicates  a  point  on  the  same  streamline  where 
the  pressure  is  computed. 

If  a  propeller  is  operating  in  a  uniform  onset  flow  with  only  an  axial  compo¬ 
nent  and  with  no  other  sources  of  disturbance,  i.e.,  the  flow  condition  for  all 
the  experimental  measurements  correlated  in  this  report,  the  pressure  will  be: 

p  =  -  yp{V-V  -  V2  -  r2^2}  +  Poo  (10) 

where  V  =  V^i  +  rfie^  +  V^,  and  p^  is  the  pressure  at  any  point  far  upstream  of  a 
propeller. 

We  define  the  pressure  coefficient  as 

c  -  -  -  4  (V-V  -  V2  -  A2)  (11) 

P  *>VR  VR 


where 


V  = 


V  = 
— w 

V  = 
~P 


V  = 
— o 


The  subscript,  A, 


where  is  a  reference  speed.  In  PSP,  three  options  are  given  for  V^;  one  is  the 
local  inflow  speed  to  the  blade  section,  /y2  +  (2mnr)^  »  the  other  two  options  are 

the  local  inflow  speed  at  r  =  0.7R  and  the  ship  speed. 


RESULTS  AND  DISCUSSIONS 

The  procedure  to  calculate  the  pressure  distribution  presented  In  the. pre¬ 
ceding  section  has  been  applied  to  the  following  configurations: 

1.  Two-dimensional  airfoil  sections;  flat  plate,  NACA  0012  and  NACA  4412 
sections 

2.  IHI  Propeller  MP  282 

3.  DTNSRDC  Propeller  4718 

4.  NSMB  Model  Propeller 

5.  DTNSRDC  Propeller  4118 

6.  DTNSRDC  Propeller  4498 

The  geometric  characteristics  of  these  propellers  are  summarized  in  Table  1.  The 
predictions  by  PSP  are  compared  with  either  experimental  measurements  or  other 
theoretical  predictions. 

TWO-DIMENSIONAL  AIRFOIL  SECTIONS 

In  order  to  test  the  validity  of  the  discrete  vortex/source  lattice  method  for 
pressure  computation,  a  computer  program,  F0IL2D,  was  developed  for  computing  the 
pressure  distribution  on  two-dimensional  airfoil  sections.  F0IL2D  has  all  the 
ingredients  of  the  discrete  vortex/source  lattice  method  except  for  three-dimen¬ 
sional  effects. 

Figure  7  shows  the  comparison  of  the  predicted  pressure  distribution  on  a  flat 

plate  with  an  angle  of  attack  a  =  4  degrees  bv  F0IL2D  with  the  analytical  solution  in 

Reference  20.  In  Figure  8,  the  pressure  distribution  is  compared  for  the  analytical 
20 

solution  and  the  F0IL2D  predictions  on  the  NACA  0012  section  at  zero  angle  of 
attack.  In  Figure  9,  measured  pressure  distribution  on  NACA  4412  section  at 
a  *  6.4  degrees  is  compared  with  predictions  by  different  methods  including  F0IL2D. 
Agreements  between  the  predictions  by  F0IL2D  and  experiments  as  well  as  those 
between  the  former  and  other  prediction  methods  are  excellent  for  two-dimensional 
shapes . 

IHI  MODEL  PROPELLER  MP  282 

The  open-water  performance  and  the  pressure  distribution  were  computed  on  the 
Ishikawajima-Harima  Heavy  Industries  (IHI)  large  model  propeller  MP  282  operating 
in  uniform  flow.  The  diameter  of  the  propeller  is  0.95  m.  This  propeller  has 
radially  varying  meanline  and  thickness  distribution.  The  predictions  are 


compared  with  the  experimental  measurements  performed  using  individual  tubes  to 
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a  hub  sensor  made  at  IHI  Ship  Model  Basin.  * 


The  open-water  performance  was  calculated  and  compared  with  experimental 
results  in  Figure  10.  The  predicted  K^,  values  are  in  excellent  agreement  with 


experimental  measurements.  The  predicted  values  are  about  5  percent  greater 


than  the  experimental  values  over  the  range  of  advance  coefficients. 

The  pressure  distribution  on  the  blade  of  Propeller  MP  282  was  calculated  for 
two  different  J  values;  J  =  1.054  and  J  =  1.163.  The  pressure  coefficients  were 
calculated  on  both  the  suction  and  the  pressure  sides  at  selected  radii  (r/R  =  0.6 
0.7,  0.8,  0.9). 


In  Figures  11  and  12,  the  experimental  measurements  and  the  predicted  are 


compared  at  J  =  1.054  and  J  =  1.163,  respectively.  The  experimental  measurements 
were  made  at  a  Reynolds  number,  R  =  1.9  x  10^.  The  calculated  pressure  coeffi¬ 


cients  are  in  good  agreement  with  measurements  on  the  pressure  side  except  near 
the  leading  edge,  but  generally  overpredict  the  suction  side  pressure.  The  agree¬ 
ment  at  the  reduced  J  value  is  better  than  that  at  the  increased  J  value.  In 
general,  the  predicted  values  are  in  satisfactory  agreement  with  the  experimental 
measurements  throughout  the  radius  at  the  two  different  J  values. 


In  Figure  13,  the  oil-film  test  results  reconstructed  from  the  photographs  in 


Reference  21  are  shown  at  two  Reynolds  numbers;  1.1  x  10^  and  2.6  x  10^.  The  oil 


film  illustrates  the  surface  streamlines  on  both  sides  of  the  blade.  At  the 


reduced  R^  condition,  the  flow  patterns  on  the  suction  side  have  significantly 
reduced  shear  stress  over  the  forward  part  of  the  blade  and  a  clear  separation 
occurs  slightly  past  midchord.  On  the  pressure  side,  reduced  shear  regions  occur 
toward  the  leading  edge  and  some  indication  of  a  leading-edge  laminar  separation 
bubble  occurs  at  both  Reynolds  numbers. 

No  surface  flow  patterns  are  presented  in  Reference  21  for  the  test  R^  of 
1.9  x  10^.  However,  judging  from  the  measured  pressure  coefficients  shown  in 
Figures  11  and  12,  it  is  possible  that  separation  occurred  near  0.7  fraction  of 
chord  on  the  suction  side  and  at  the  leading  edge  on  the  pressure  side  in  the  form 
of  a  bubble.  Such  separation  would  explain  the  suction  peak  on  the  pressure  side 
near  the  leading  edge  and  the  pressure  peak  measured  at  0.7  radius  at  0.7  fraction 
of  chord  (measurements  were  not  made  at  a  similar  chordwise  position  at  other 
radii).  It  is  further  hypothesized  thai.  the  suction  side  separation  is  a  thin 


layer  with  only  minor  influence  on  the  pressure  away  from  the  separation  line. 
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Previous  data  for  this  propeller  indicated  that  the  pressure  at  the  point  in 
question  exhibited  the  same  property  as  a  function  of  Reynolds  number. 

DTNSRDC  PROPELLER  4718 

The  steady  pressure  distribution  was  calculated  on  the  surface  of  DTNSRDC 
controllable-pitch  Propeller  4718  at  the  design  advance  coefficient,  J  =  0.75. 

The  propeller  has  three  blades  with  diameter  of  2  feet  (0.61  meters),  EAR  of  0.44, 
and  tip  skew  of  20  degrees.  In  Figure  14,  the  predicted  pressure  coefficients  are 
compared  with  experimental  measurements"*  on  the  blade  surface  with  the  propeller 
operating  in  uniform  flow  at  three  different  radii;  r/R  =  0.5,  0.7,  and  0.9.  The 
pressure  was  measured  by  transducers  mounted  on  both  sides  of  the  blade  surface. 

The  experimental  values  were  measured  at  six  different  R  values  ranging  from 

6  6  ^ 

2.5  x  10  to  4.63  x  10  at  the  design  J.  For  this  range  of  Rn>  the  flow  on  the 

surface  should  be  fully  turbulent  so  that  the  pressure  distribution  would  be  nearly 

independent  of  R^ .  However,  the  experimental  measurements  showed  substantial 

variations  for  different  R^  values.  The  variation  is  more  pronounced  on  the 

suction  side  than  on  the  pressure  side.  In  general,  the  pressure  coefficients 

increase  with  increasing  R^  values.  The  measured  pressure  coefficients  shown  in 

Figure  14  represent  the  average  values  over  the  range  of  R^.  Unpublished  flow 

visualization*  of  the  surface  streamlines  showed  no  anomolous  flow  over  the  blades. 

At  r/R  =  0.7,  the  computed  values  are  in  reasonable  agreement  with  experi¬ 
mental  values.  However,  the  agreement  at  the  other  two  radii  is  not  as  good  as 
that  at  r/R  =  0.7.  At  r/R  =  0.5,  the  experimental  results  show  some  irregular 
peaks  at  x£  =  0.12  and  0.5  on  both  sides.  Jessup"*  explained  that  some  of  these 
irregularities  in  the  measured  values  might  be  partially  attributed  to  the  effect 
of  the  relatively  large  fairwater  and  hub. 

More  recently,  Jessup*  measured  blade  pressure  on  the  same  Propeller  4718 
using  another  technique.  In  this  experiment,  Jessup  measured  the  pressure  distri¬ 
bution  only  on  the  suction  side  at  two  radii,  r/R  =  0.5  and  0.8  at  the  design 
J  =  0.751.  In  Figure  15,  the  two  sets  of  experimental  values  are  compared  with  PSP 


*Private  communication  from  S.  Jessup,  DTNSRDC,  Code  1544. 
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predictions.  Although  both  experimental  measurements  showed  Reynolds  number 
effects,  the  correlation  between  PSP  predictions  and  the  new  experimental  values 
improved  substantially. 

The  larger  discrepancy  at  r/R  =  0.9  on  the  suction  side  may  be  due  to  real 
flow  effects.  However,  experimental  inaccuracy  demonstrated  at  r/R  =  0.5,  or  the 
coarse  modeling  for  the  global  solution  in  the  analytic  treatment  of  the  flow  in 
that  region  can  also  be  a  possible  source  of  the  discrepancy. 

NSMB  MODEL  PROPELLER 

The  steady  pressure  distribution  was  computed  on  the  NSMB  model  propeller  at 
J  =  0.4  and  0.6  to  correlate  with  experimental  data.^  This  propeller  has  simple 
geometric  characteristics;  no  rake,  no  skew,  and  a  single  section  shape  over  the 
radius.  It  was  originally  designed  for  bubble  cavitation  investigation. 

The  open-water  performance  was  calculated  and  compared  with  experimental  re¬ 
sults  in  Figure  16.  The  predicted  K^,  and  are  in  excellent  agreement  with  the 
experimental  values  for  the  range  of  J  values  except  for  very  reduced  ones. 

In  Figures  17  and  18,  the  predicted  pressure  coefficients  are  compared  with 
experimental  measurements  obtained  at  NSMB ^  at  J  *  0.4  and  0.6,  respectively,  at 
five  different  radii;  r/R  =  0.4,  0.5,  0.6,  0.7,  and  0.8.  At  r/R  =  0.9,  only  pre¬ 
dicted  values  are  plotted  since  the  pressure  was  not  measured  at  that  radius.  The 
pressure  was  measured  by  transducers  mounted  within  both  sides  of  the  blade  surface. 

The  experimental  measurements  were  made  twice  within  a  six-month  period  in 
order  to  assess^  the  repeatability.  The  two  series  of  experimental  results  are 
shown  in  Figures  17  and  18.  The  repeatability  is  generally  good.  The  predicted 
values  on  both  sides  are  in  good  agreement  with  experimental  measurements  at  both 
J  values. 

COMPARISON  WITH  OTHER  PROCEDURES 

In  1978,  the  ITTC  Propeller  Committee  surveyed  existing  prediction  methods  for 
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pressure  distribution  on  the  propeller  blade  surface.  They  compared  the  pre¬ 
dictions  made  by  various  methods  from  sixteen  participating  institutions  throughout 
the  world.  The  propeller  selected  for  the  comparative  calculations  was  DTNSRDC 
Propeller  4118,  a  three-bladed  research  propeller  tested  thoroughly  at  DTNSRDC  for 
open-water  performance,  cavitation,  and  unsteady  forces. 


In  Figure  19,  the  predictions  by  PSP  are  compared  with  other  predictions 
represented  by  the  envelope  covering  all  the  predicted  results  at  the  design 
advance  coefficient,  J  =  0.833.  The  predicted  values  by  PSP  are  within  the 
envelope  of  the  predictions  by  other  methods. 

In  Figure  20,  comparisons  are  made  between  predictions  by  PSP  and  by  a  lifting 
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surface  method  presented  by  Brockett  for  a  propeller  similar  to  DTNSRDC  Propeller 
4498  at  J  =  0.888.  The  propeller  is  warped  with  72  degrees  warp  angle  at  the  tip. 
The  section  meanline  is  similar  to  the  NACA  a  =  0.8  meanline. 

The  predictions  made  by  the  two  different  methods  are  in  good  agreement  at 
r/R  =  0.254,  but  the  discrepancies  increase  toward  the  tip  region,  as  it  did  for 
the  experimental  data  of  Jessup. 

SUMMARY  AND  CONCLUSIONS 

The  discrete  vortex/source  lattice  lifting  surface  method  has  been  used  for 
the  prediction  of  steady  pressure  distribution  on  a  rotating  propeller  blade  sur¬ 
face.  A  computer  code,  PSP,  has  been  developed  by  extending  the  existing  propeller 
global  analysis  program,  PSF,  developed  at  M.I.T. 

For  pressure  computations  on  the  propeller  blades,  the  velocity  jump  across 
the  vortex/source  sheet  must  be  carefully  treated  and  include  the  effects  of  both 
the  spanwise  and  chordwise  vortices.  In  PSP,  the  effect  of  the  chordwise  vortices 
at  the  pressure  point,  the  midpoint  of  each  spanwise  vortex,  was  accounted  for  by 
interpolating  from  the  four  adjacent  chordwise  vortices. 

Comparisons  of  the  predictions  by  PSP  with  experimental  measurements  ana  pre¬ 
dictions  by  other  methods  on  selected  model  propellers  generally  showed  good 
correlations.  The  correlations  near  the  tip  region,  especially  for  skewed  pro¬ 
pellers,  i.e..  Propellers  4718  (20  degrees  tip  skew)  and  4498  (72  degrees  tip  skew) 
are  not  as  good  as  those  for  the  inner  region.  Possible  explanations  may  be  that 
near  the  tip  region  of  skewed  propellers,  viscous  effects  may  be  large  or  that  the 
current  numerical  modeling  in  lifting  surface  representations  may  not  be  accurate 
enough . 

*The  predictions  by  Brockett  shown  in  Figure  20  are  taken  from  Figure  A 
(linear  3D  method)  in  "Discussions  and  Authors'  Closures"  section  of  Reference  9. 


RECOMMENDATIONS 


Based  on  the  investigations  made  in  the  present  work,  the  following  studies 
are  recommended  in  order  to  further  improve  the  current  prediction  method: 

1.  The  improved-accuracy,  tip-flow  part  of  the  PSF  should  be  used  for  the 
prediction  of  the  pressure  distribution  near  the  tip  region.  The  tip  flow  is  very 
complicated  and  of  practical  importance,  and  yet  the  prediction  near  the  tip  region 
is  not  as  good  as  that  for  the  inner  region.  Since  the  tip  flow  model  contains 

a  finer  lattice  arrangement  than  does  the  global  flow  model,  the  tip  flow  solution 
is  expected  to  give  more  accurate  results.  The  modification  to  the  tip  flow  part 
for  pressure  calculations  is  straightforward. 

2.  In  order  to  be  able  to  predict  viscous  phenomena  such  as  suction-side 

separation  or  leading-edge  laminar  bubble  separation  that  is  frequently  observed 

in  experiments  with  model  propellers,  suitable  analytical  and  numerical  analysis 
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should  be  undertaken.  Some  initial  efforts  in  this  area  have  been  undertaken. 

3.  In  order  to  further  assess  the  validity  and  limitation  of  the  current 
procedure,  comparative  calculations  are  recommended  with  other  theories  for  a 
wide  range  of  propellers  and  operating  conditions. 

4.  Parametric  calculations  of  propeller  characteristics  of  practical 
importance  such  as  cavitation  inception,  boundary  layer  development,  and  blade 
stress  should  be  undertaken. 
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APPENDIX 

PRESSURE  IN  A  MOVING  FRAME  OF  REFERENCE 


and 


In  a  moving  frame  of  reference,  Euler's  equation  of  motion  of  an  inviscid 
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incompressible  fluid  can  be  expressed  as  follows  (see  Batchelor  ): 


DV 

Dt 


-  —  Vp 

P 


+  F 


(  d^r  dfl 
dt 


dt 


xr  +  2flxV  +  flx  (fix  r)' 


(12) 


D  9 

material  derivative  defined  by  —— ■  =  +  V  •  V 

Dt  d  t 

total  velocity  with  respect  to  the  moving  reference  frame 
fluid  density 
pressure 

body  force  per  unit  mass 

position  vector  of  the  origin  of  the  moving  frame 

angular  velocity  of  the  moving  frame  about  the  origin 
position  vector  of  a  field  point  in  the  moving  frame 

The  last  two  terms,  2S1  x  V  and  fl  x  (S)  x  r) ,  are  called  the  Coriolis  force  and  the 

centrifugal  force,  respectively. 

If  we  take  r  =0  and  Q  =  constant.  Equation  (12)  becomes: 

— o  — 

3V  j 

+  (V*V)V  =  -—  Vp  +  F  -  2fi  x  V  -  SI  x  (fl  x  r)  (13) 

dt  P 

Assuming  that  the  gravity  force  is  the  only  body  force  acting  on  the  fluid,  one 
can  express  F  by: 

F  =  V (-gy  )  (14) 

—  o 

where  g  is  the  gravitational  acceleration  and  yQ  is  the  vertical  coordinate  in  the 
nonrotating  coordinate  system  (x^.y^z^)  as  shQVm  Figure  1.  It  is  to  be  noted 
that  this  term  is  time-dependent  in  the  rotating  frame  of  reference. 

Utilizing  Equation  (14)  and  the  vector  identities, 


where 


D 

Dt 


Y  = 

p  = 

p  = 
F  = 


r 

—o 


n  = 

r  = 
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yV(V*V)  =  (V*V)V  +  V  x  (V  x  V) 


-Q  x  (Q  x  r)  =  -fvjft  x  r|^, 


one  can  express  Equation  (13)  as  follows: 

V  j-f  V«V  +  -jj-  -  y  |ft  x  r|2  +  gyQ 


+  V  x  (to  +  29.) 


Here,  co  xs  the  vorticity  in  the  fluid  measured  in  the  moving  frame  of  reference. 

Now,  consider  a  propeller  rotating  at  a  constant  angular  velocity,  fi  =  -fti_ 

(a  right-hand  rotation  propeller,  see  Figure  1),  in  an  axisymmetric  wake  of  a  ship, 
where  the  flow  is  steady  in  the  rotating  frame  of  reference.  The  ship  wake  velo¬ 
city,  V^,  can  be  expressed  in  polar  components  as: 


V  =  V  i  +  V  e  +  Vnen  (18) 

— w  x—  r— r  0-6 

where  V^,  and  are  radially  varying  axial,  radial,  and  tangential  components, 
respectively.  It  is  assumed  that  the  variation  of  the  ship  wake  velocity  in  the 
radial  direction  is  small. 

In  addition  to  the  ship  wake  velocity  we  assume  that  there  exists  another 
axisymmetric  disturbance  velocity  component,  V^,  that  is  introduced  locally  by 
nearby  appendages  or  other  lifting  surfaces: 

V  =  V  i+V  e  +  V  .e_  (19) 

— o  ox—  or—  r  00— 0 

Then  the  total  velocity,  V,  in  a  cylindrical  coordinate  system  rotating  with 
the  propeller  can  be  expressed  as  follows: 

V  =  V  +  rfte-  +  V  +  V  (20) 

—  -w  ^0  — p  — o 

where  V  is  the  perturbation  velocity  due  to  the  presence  of  the  propeller. 

'P 
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In  the  rotating  coordinate  system,  the  vorticity  vector,  0),  can  be  expressed 
as  the  sum  of  two  terms: 

Cl)  =  CD  +  CD  (21) 

—  — r  — w 

where  cor  is  the  vorticity  due  to  the  rotation  of  the  coordinate  system  and  c^  is 
the  vorticity  in  the  inertial  reference  frame.  From  the  definition  of  vorticity, 
it  can  be  shown  that: 

(D  =  V  x  (rSiO  =  -20  (22) 

r  — 0  — 

If  we  let  _r  =  xi  +  re^(8),  |fl  x  r|  term  on  the  left-hand  side  of  Equation  (17) 
will  be: 

|Q  x  r|  =  rfl  (23) 

Substituting  Equations  (21)  to  (23)  into  Equation  (17)  with  the  assumption 
of  the  steady  flow,  we  have: 

V  Jy  V*V  +  -2-  -  y  (rQ)2  +  gyQj  =  V  x  c^  (24) 


By  integrating  Equation  (24)  along  a  path  in  the  flow  between  two  arbitrary 
points,  A  and  B,  we  obtain  the  Bernoulli  equation: 


hb  -  HA  " 


(V  x  0)  )  •  dr 
—  -w  — 


where  H(r, t)  is  sometimes  called  the  Bernoulli  head  and  is  defined  by: 


H(r, t)  =  -L  vv  +  -2-  -  -f  (rii)2  +  gy 
—  l - p  z  o 


If  we  take  the  integral  path  dir  along  a  streamline  or  a  vortex  line,  i.e., 
parallel  to  V  or  o^,  respectively,  the  integral  in  Equation  (25)  vanishes  since  the 
dot  product  in  Equation  (25)  is  equal  to  zero.  It  then  follows  that  the  Bernoulli 
head  is  constant  along  a  streamline  or  a  vortex  line. 

By  taking  a  reference  point,  A,  as  a  point  along  the  streamline  far  upstream 
of  the  propeller  where  the  propeller  perturbation  velocity,  V  ,  and  the  other 
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disturbance  velocity,  V  ,  are  negligibly  small,  the  Bernoulli  constant,  H.,  will 

O  A 

be: 

HA  ■  [t  v* +  f  +  «\,]A 

Then  the  pressure  at  an  arbitrary  point  in  the  fluid  can  be  expressed  as: 

Pa'TP  {-*-  "  (Va  “  '  P8(yo  ‘  yoA}  +  PA  (28) 

where  V  =  +  V  +  Vq,  and  the  subscript  A  indicates  a  point  on  the  same 

streamline  (or  vortex  line)  where  the  pressure  is  computed.  The  effect  of  gravity 
-pg(y  -  y  )  in  Equation  (28)  gives  rise  to  a  once-per-revolution  periodic  varia- 

O  O  A 

tion  in  the  pressure  in  the  rotating  coordinate  system.  Since  this  term  does  not 
contribute  to  the  mean  pressure  and  the  loading,  it  is  not  considered  in  the  pre¬ 
sent  study.  However,  this  term  may  be  important  when  cavitation  inception  is  of 
interest. 

For  a  uniform  onset  flow  (potential  flow)  with  only  an  axial  component  and 
with  no  other  disturbance  than  the  propeller  itself,  i.e.,  the  flow  condition 
applicable  to  all  the  experimental  measurements  correlated  in  this  report,  the 
pressure  equation  becomes  even  simpler: 

p  .  -  |p|v.V  -  V2  -  r2n2^+  (29) 

where  V  =  V^i  +  rfle^  +  V^,  and  p^  is  the  known  pressure  far  upstream.  In  this 
case,  the  Bernoulli  head  is  constant  everywhere  in  the  fluid  since  there  is  no 
vorticity  in  the  flow  (see  Equation  (24)). 


Figure  A  -  Trailing  Wake  Geometry  after  Wake  Alignment 
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6  -  Schematic  Representation  of  the  Effect  of  Chordwise 
Vortices  on  the  i^  Pressure  Point 
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ure  9  -  Pressure  Distribution  on  NACA  4412  Section 
at  a  =  6.4  Degrees 
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Figure  10  -  Open-Water  Performance  of  IHT  Model 
Prooeller  MP  282 
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Figure  13  -  Surface  Flow  Patterns  by  Oil-Film  Test  on  T H I 
Model  Propeller  MP  282 
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Figure  14  -  Pressure  Distribution  on  DTNSRDC  Propeller 

4718  at  J  -  0.75 


DTNSRDC  PROPELLER  4718,  J  =  0.75 


Figure  15  -  Pressure  Distribution  on  Suction  Side  of  DTNSRDC  Propeller  4718  at 
r/R  =  0.5:  Comparison  of  Prediction  and  Two  Sets  of  Experiments 
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Figure  16  -  Open-Water  Performance  of  NSMB  Model  Propeller 
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Figure  18  (Continued) 
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Figure  19  -  Pressure  Distribution  on  DTNSRDC  Propeller  4118  at 
J  =  0.833:  Comparison  with  ITTC  Propeller  Committee  Survey 
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